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Abstract: 

We study the three dimensional Navier-Stokes equation with a random Gaussian force 
acting on large wavelengths. Our work has been inspired by Polyakov's analysis of steady 
states of two dimensional turbulence. We investigate the time evolution of the probability 
law of the velocity potential. Assuming that this probability law is initially defined by a 
statistical field theory in the basin of attraction of a renormalisation fixed point, we show 
that its time evolution is obtained by averaging over small scale features of the velocity 
potential. The probability law of the velocity potential converges to the fixed point in the 
long time regime. At the fixed point, the scaling dimension of the velocity potential is 
determined to be — |. We give conditions for the existence of such a fixed point of the 
renormalisation group describing the long time behaviour of the velocity potential. At 
this fixed point, the energy spectrum of three dimensional turbulence coincides with a 
Kolmogorov spectrum. 
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Turbulence is one of the crucial puzzles of theoretical Physics. Despite its practical 
relevance, there is still no real thorough understanding of the principles of turbulence. A 
few years ago Polyakov suggested using conformal field theories to study turbulence in 
two dimensions . Polyakov's idea is that the scale invariant regime of turbulence can 
be described by a conformal field theory. In three dimensions, perturbative renormalisa- 
tion transformations in momentum space have been used to derive an infrared fixed point 
describing the long distance physics of the Navicr-Stokes equationt^'*^'"^]. Similarly, prob- 
abilistic solutions of the Navier-Stokes equation with no forcing term have been provided 
by self-similar probability laws for the velocity field t^'^L We shall use properties of the 
Navier-Stokes equation under scaling transformations in order to study the time evolution 
of solutions. In particular, our analysis will combine non-perturbative statements about 
the renormalisation group and results about 3d conformal field theories applied to the 
Navier-Stokes equation. 

The motion of a fiuid is described by the macroscopic velocity field. We shall suppose 
that the macroscopic velocity field Va is regularised at small distances by a cut-off a. This 
cut-off is related to the viscosity by 



- (1) 

where r is a time characterising the energy decay when a stirring force is absent. For 
incompressible fluids, the velocity field is divergenceless: 

^.Va = 0. (2) 

This implies that the velocity field can be written 

Va = curl ijja (3) 
where ipa is the velocity potential. We will be interested in the vorticity 

Ua = curl Va- (4) 

The vorticity satisfies the vorticity Navier-Stokes equation 



dt 



+ {Va.V)uJa = {cOa.V)Va + l^AcJa + Ma. (5) 



where M — curlF is the vorticity forcing term and F is the velocity forcing term. In order 
to discard boundary effects, we shall only be interested in the case where the fluid fills up 
all space. Furthermore, we shall require the vorticity to decrease at infinity. The initial 
conditions and the external stirring are chosen to be regularised random fields where the 
ultra violet cut-off is given by the scale a. In particular, the velocity forcing term will 
be supposed to be Gaussian with a white noise temporal dependence and an action over 
large distances. Taking into account the incompressibility condition, its 2-point correlation 
function is 

< FaAk,uj)FaAk',uj') >= WoSL{k){S^f3-^)S{k + k')S{u; + uj') (6) 
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in Fourier space. The constant Wq determines the amphtude of the forcing while its k- 
dependence is defined by SL{k). We suppose that Sl is a positive smoothing function of 
width L~^, i.e. SL{k) = L^S{Lk) where Sl has width and is centered around zero, it 
is also normalised by J d^kS{k) — 1. The width is supposed to be arbitrarily small. Notice 
that the spectrum'* of the stirring force F is k^SL{k) which is peaked around k L~^. As 
a is the only intrinsic scale, the scale L is related to a by 

L = Ka (7) 

where k is a large number depending on the forcing term. For a given initial configu- 
ration and a given realisation of the stirring force, solutions of (5) will solve the initial 
value Cauchy problem. Starting from an initial probability law dVo{ip), one can solve (5) 
for a given realisation of the stirring force. This defines the conditional probability law 
dV{ipa{x, t)\Fa{x, t)) expressing the probability law of ipa{x, t) for a given realisation of the 
stirring force Fa- This probability law is such that the Navier-Stokes equation is valid in 
all correlation functions 

d / X / \ 

— < UJa{Xi,t)...UJa{Xn,t) > = 

n (8) 
22 < '^a{xi, t)...{-{Va.V)uJa + (Wa-V)t'a + I^AuJa + Ma) ...UJa{Xn, t) > 

for a given realisation of the stirring force Fa. These equations are the Hopf equations. The 
conditional probability law dV{il)a{x, t)) is random when seen as a functional of the 

stirring force. After averaging over the forcing term realisations, one gets the probability 
law of the velocity potential at time t 

dViM^^t)) = J dV{Mx,t)\Fa{x,t))®dfx{Fa{x,t)) (9) 

where diJ,{Fa{x, t)) is the Gaussian probability law of the stirring force at time t. The main 
object of the present work is to characterise the long time behaviour of dV{'ipa)- 

We shall suppose that the conditional probability at time t is defined by a statistical 
field theory 

dV{Mx,t)\Fa{x,t)) = ^di^exp-St,Ma,Fa) (10) 

where the partition function Z is a normalisation factor, dif) stands for the functional inte- 
gral symbol and St^rii^i - Fa) is an effective action at time t depending on the regularised 
fields. By dimensional analysis, the action is a homogeneous function of ^. After integrat- 
ing over the stirring force, one obtains the effective action describing the probability law 
of the velocity potential at time t 

dViipa) = ^#exp-S't,^(V'a) 

(11) 

e^p-St,r{'ipa) = / eW-St,r{'ipa,Fa)dl^{Fa{x,t)) 



Practically, S(k) can be taken to be TT 2 exp— A;^, in that case k'^SL{k) has a maximum at 
k = L-^. 
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We shall suppose that the velocity potential is a field with scaling dimension d^p. This 
dimension is for instance well-defined if the effective field theory (11) is the basin of at- 
traction of a fixed point of the renormalisation group. In that case the scaling dimension 
of ip becomes its conformal dimension. In this communication, we shall show that the 
time evolution of solution of (9) starting from the initial conditions (11) is determined by 
renormalisation transformations. The evolution of solutions of (9) is given by 

/(12) 
dijex.p-St,x-T^(^x-^aix,t)). 

This is the precise statement that the time evolution of solutions follows renormalisation 
trajectories. The first equation implies that the velocity potential lies on the renormalisa- 
tion trajectory of a field of dimension d^, i.e. the solution i/jaiXx, X'^t) of the Navier-Stokes 
equation at time A^t is simply the renormalised field tpaiXx^t) of ipx-^aix-it)- The second 
equation entails that the effective action at time X^t is obtained after integrating over the 
fluctuations in the range [A~-^a, a]. This renormalisation transformation is performed in 
momentum space integrating over modes in [a~^ . Xa~^]. The time rescaling exponent T 
and the scaling dimension of if) are uniquely determined 

4 (13) 
d^ = -- 

Notice that the scaling dimension of a; is = | guaranteeing that the vorticity decreases 
at infinity. 

Let us now prove these results. They depend crucially on properties of the stirring 
force Fa- Notice that the spectrum of the stirring force SL{k) has length dimension —3. 
Moreover, one obtains 

SLij) = X'Sx-.Lik) (14) 
Using (6), one can deduce the following scaling relation 

< Faij, ^),Fa{j, ^) >= A«+^ < Fx-.a{k,u),Fx-.a{k',u;') > (15) 
As the stirring force is Gaussian, this relation between 2-point functions yields an equality 

Fa{Xx,X^t) = X=^Fx-ia{x,t), (16) 

i.e. these two random variables have the same probability law. In other words, one has 
the equality 

J f{X^Fa{Xx,X'^t))df,{Fa{Xx,X'^t)) = J f{Fx-ia{x,t))df,{Fx-ia{x,t)) (17) 
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for any functional /. This fact will allow us to relate the probability law of the stream 
function at two different times when taking the average of (19). Let us rewrite the Navier- 
Stokes equation after rescaling a —>■ X~^a and r —>■ X~^t 



dt 



At time t, the initial condition (10) defined by St,^{il^\-i{x^t),Fx-ia{'^it)) is a solution 
of (18) where the range of the integration over modes for any correlation function is now 
[0, Aa"^]. Using (12) and (16) one can substitute A^^V'a(Ax, X^t) and A^+^Ma(Ax, X^t) for 
'4^\-'^a{x, t) and Mx-ia{x, t) in each correlation function. Notice that the fields in each cor- 
relation function have now modes only in the range [0, a~^] . The integration over the modes 
in the range [a~^ ^ Xa~^] only affects the Boltzmann weight exp —Si \-T^{'^x-ig_{x, t), Fx-ig^{x, t)).| 
Performing the integration of the Boltzmann weight over these modes gives the renor- 
malised Boltzmann weight 

exp -SxTt,r{^a{Xx, X^t), X^ Fa{X^t, Xx)) = / di; exp -St,^{'ipx-^a{x,t), Fx-ia{x,t)) 

(19) 

The probability law of the velocity potential (12) is obtained after averaging over the 
Gaussian stirring force and using (17). The resulting equations for the correlation functions 
only involves the fields V'a(Aa;, X'^t) and Ma(Ax, X^t). These equations correspond to 

^""^ + X^^''*-^{Va.dx.)Ua = X'+'^-^icOa.dxa.H + udlcOa + X''-^"'^ Ma (20) 



d{xn) 

when inserted in correlation functions. The derivatives are taken with respect to X^t and 
Xx. We now require that (20) coincides with (5) at time X^t and coordinates Xx. This is 
achieved if 

T 

d^ = —1 , , 

^ 2 (21) 

d^ = T — 2. 

One can then deduce (13). This proves that (12) and (13) specify the time evolution of 
the probability law of the velocity potential. 

Starting from a field theory in the basin of attraction of a fixed point of the renormal- 
isation group, this probability follows renormalisation trajectories. It therefore converges 
to the probability law specified by the fixed point. The scaling dimension — | of the ve- 
locity potential becomes its conformal dimension at the fixed point ( fixed points of the 
renormalisation group are conformal theoriest^'). We shall now derive further conditions 
for the existence of such a fixed point. As the long time regime is linked to the small j 
region by (12), the fixed points satisfy (18) in the limit when A goes to infinity. Using (13), 
one can see that the influence of the viscosity becomes negligible in the long time regime 
as T — 2 = — I . Similarly, the stirring force can be easily dealt with. Indeed, notice that 

S.{k) = {^rS{^) (22) 
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converges to zero as A — > oo. This imphes that the 2-point correlation of F^-i^ goes to zero 
and therefore the Gaussian stirring force vanishes in the long time regime. Assuming that 
the fixed points represent steady states, the time derivatives of each correlation function 
vanish as well. We therefore find that fixed points are characterised by the vanishing of the 
non-linear terms in the limit A ^ oo. The non- linear terms can be readily evaluated when 
the cut-off is rescaled to zero using properties of short distance expansions. Products of 
fields become singular when the cut-off is removed. For fields called quasi-primary fields 
the result of the product can be expanded in a power series in ^. A vanishing limit is 
obtained if the leading term has a positive exponent. We shall suppose that the velocity 
potential becomes a quasi-primary field at the fixed point. In that case, the nonlinear 
terms read 

ea/37^7(K-ia-V)t;A-ia)/3 ~ (^)'^-'-''* V'2,a,A-ia + - (23) 

where ^2 is the conformal dimension of the leading pseuso- vector 7/^2, a in the expansion of 
the non-linear terms. We can therefore conclude that the fixed points are characterised by 
the inequality 

d2 > 4 + 2d^ (24) 

This generalises a similar inequality obtained by Polyakov in two dimensions. We can 
now state our main result. Conformal field theories satisfying (24) such that the velocity 
potential is a quasi-primary field of dimension — | describe the long time regime of 3d 
turbulence. 

As a consequence, notice that -0 cannot be identified with ■02 as (24) is not satisfied 
in that case. This implies that 

< V'mV'^^V'p >= 0, (25) 
the three point function of the velocity potential vanishes identically at the fixed point. 
It is extremely difficult to construct explicitly three dimensional conformal field theories 
describing an infrared fixed point of the renormalisation group. We shall suppose that 
such theories exist and deduce consequences on the energy spectrum of the solutions of the 
Navier-Stokes equation. In particular, we can calculate the 2-point function of the velocity 
potential and then the energy spectrum in the long time regime. The energy spectrum is 
easily related to the Fourier transform of the two point function 

Ea{k) = A-Kk" J d^X < Va{x)Va{0) > exp —2'Kik.x. (26) 

In the long time regime, the energy spectrum behaves as 

Ea{k)^k-i. (27) 

This spectrum is valid for A; <^ ^. The exponent is given by 2d^ + 1 = — | The spectrum 
coincides with a Kolmogorov spectrumt^l 

^ Our results are compatible with the perturbative calculations of Refs. [2,3,4] . In this case, the stirring 
force is specified by < (/c, cu)F„,^(/c', w') >- ^\ y Xa-^{k){5^p-^)6{k+k')5{uJ+uj') 

where X^-i (/c) is equal to one for k < 0,^^ and zero otherwise. The mass term m — is small but 

guarantees that there are no infrared divergences when k — 0. The same analysis as in the text shows 
that di[, = — and the spectrum behaves like k s in the range k < for < y < 5. 
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The solutions of the Navier-Stokes equation (12) have some sahent features. Let us 
notice that the evolution of the probability law of the velocity potential (12) makes explicit 
the breaking of time reversal invariance already present in the Navier-Stokes equation. 
Indeed, the probability law at time X^t is obtained from the probability law at time 
t by averaging over the small scale features of the velocity potential. These small scale 
features are therefore lost in the process of evolution. In particular, knowing the probability 
law at time X^t does not enable one to deduce the probability law at time t. Another 
consequence is that the renormalisation trajectories converge to a non-unitary conformal 
field theory as c?^ < 0. This result has also been obtained in 2dt^^. It reinforces Polyakov's 
idea that turbulence is a flirx state and not an equilibrium state described by an unitary 
theory. Moreover, at the fixed point the energy spectrum is determined to be a Kolmogorov 
spectrum. This result highly depends on the fact that the stirring force is Gaussian and acts 
on large wavelengths. Another feature of our analysis is that we have derived a spectrum 
compatible with a Kolmogorov spectrum without any cascade hypothesis t^^. Finally, let 
us note that our results can be applied in two dimensions, in particular (24) is replaced by 
Polyakov's conditiont^^^ and the energy spectrum is still given by a Kolmogorov spectrum. 

We have been able to obtain properties of solutions of the 3d Navier-Stokes equa- 
tion in the long time regime. Assuming that the probability law of the stream function 
is a statistical field theory, its time evolution corresponds to renormalisation trajectories. 
Starting from initial conditions in the basin of attraction of a fixed point under renor- 
malisation transformations, the solutions converge asymptotically to the fixed point. The 
possible fixed points are all non-unitary. In the case where the random stirring has a 
spectrum concentrated on large wavelengths, we have shown that a Kolmogorov spectrum 
is obtained in the long time regime. The existence of such a fixed point of the renormalisa- 
tion group governing the long time regime of turbulence seems physically clear in view of 
the numerous experimental as well as computational evidence in favour of a Kolmogorov 
spectrum. It is a challenging problem to construct such a field theory, 
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^ The stream function Ip satisfies ll^\-ia{x)'^\-^a{x) ~ ij)'^^ ^'^'^ '^2,\-^a{x) with ^2 > 2d^. 
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